Energy scan by <p mesons and threshold energy for the confinement-deconfinement 

phase transition 
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We argue that the ratio of <j> mesons multiplicity over cube of the mean pr is proportional to the 
degeneracy of the medium produced in ultra-relativistic heavy ion collisions. The ratio extracted 
from the existing <j> meson data in the energy range ^=6.3-200 GeV, indicate that beyond a 
threshold energy y/s th = 15.74 ± 8.10 GeV, the medium crosses over from a confined phase to a 
deconfined phase. 

PACS numbers: 25.75.-q, 25.75.Nq, 25.75.Ag 



Lattice QCD predicts [J 

that in ultra-relativistic 
heavy ion collisions, a confinement-deconfinement phase 
transition can occur, producing a new state of matter, 
Quark-Gluon Plasma (QGP). QGP is a collective state 
where color degrees of freedom become manifest over 
nuclear rather than hadronic volume. Recent experi- 
ments have produced convincing evidences for a 
confinement-deconfinement phase transition in y / s=200 
GeV Au+Au collisions at RHIC. One naturally won- 
ders, whether or not there is a threshold energy for the 
transition. One of the aims of the STAR's energy scan 
programme at RHIC is to determine the threshold en- 
ergy for the confinement-deconfinement transition 0, Q . 
STAR proposes to study nuclear collisions at y/s— 5, 7.7, 
11.5, 17.3, 27 and 39 GeV. With the existing data at 
Vi=62,130 and 200 GeV, STAR will scan a large energy 
range, ^1=5-200 GeV. 

For long, strangeness enhancement is considered as a 
signature of QGP formation In QGP environment, 
gg — > ss is abundant. If not annihilated before hadro- 
nisation, early produced strange and anti-strange quarks 
will coalesce in to form strange hadrons and compared to 
elementary pp collisions, strange particle production will 
be enhanced. However, strangeness enhancement could 
also be obtained in a purely hadronic scenario, mainly 
due to the 'volume effect' |10r4l3l |. Strangeness produc- 
tion in small volume elementary pp collisions can be 
'canonical ly' sup pressed due to 'strict' strangeness con- 
servation [U>[l3l |. In bigger volume A A collisions, lo- 
cally, strangeness conservation condition can be relaxed 
to produce strange particles. In the language of statisti- 
cal mechanics, while canonical ensemble is applicable in 
pp collisions, grand canonical ensemble is more appropri- 
ate in heavy ion collisions. Additionally, strange particle 
phase space appears to be under-saturated in elementary 
pp or peripheral heavy ion collisions [lil [l5| . 

As noted in [l6j], several unique features of (f> mesons 
make it an ideal probe to investigate medium properties 
in heavy ion collisions. They are hidden strange particle 
and are not affected by canonical suppression. Also they 
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FIG. 1: Black circles are lattice simulation 2] for s/T 3 , en- 
tropy density over cube of the temperature. The solid line is 
a parameterisation of lattice simulation of s/T 3 . 



are not affected by resonance decay, have both hadronic 
and leptonic decay channels, mass and width are not 
modified in a medium [l7j . etc. Incidentally, experimen- 
tal data for (j) meson production in nuclear collisions over 
a large energy range y / i=6-200 GeV exists. NA49 col- 
laboration measured <j> meson in 20A, 30A, 40A, 80A and 
158A GeV Pb+Pb collisions " 



171 ] . The centre of mass en- 
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ergies are, 6.3, 7.6, 8.3, 12.3 and 17.3 GeV respectively. 
STAR collaboration measured <\> mesons in Au+Au and 
Cu+Cu collisions at RHIC pB^. In Au+Au collisions 
(j) mesons are measured at three energies, y / s=62, 130 
and 200 GeV, and in Cu+Cu collisions at y/s=62 GeV 
and 200 GeV. In the following we argue that the exist- 
ing <f> meson data do indicate a threshold energy for the 
confinement-deconfinement phase transition. 

Let us note the most distinguishing feature between a 
confined and a deconfined state. Effective degrees of free- 
dom of a deconfined phase is considerably larger than 
that of a confined phase. As an example, in FigHJ a 
recent lattice simulation Q for the temperature depen- 
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FIG. 2: The black circles hydro-dynamic model simulations 
for <j> mesons mean px as a function of spatially averaged 
initial temperature < T, >. The solid line indicate that in a 
hydrodynamic model, <f> mesons mean px depend linearly on 
the average initial temperature. 



FIG. 3: The filled symbols are the ratio of <f> meson multiplic- 
ity over the cube of the mean px, normalised by participant 
number (N par t) and rapidity gap (AY), in Pb+Pb, Cu+Cu 
and Au+Au collisions, as a function of the collision energy. 
The solid line is a fit to the ratio by the analytical form for 
the step function EqO the dotted line is a power law fit. 



dence of entropy density over cube of the temperature 
(s/T 3 ) is shown. The simulation was done with two light 
flavor quarks and a heavy strange quark, with almost 
physical quark masses. s/T 3 is approximately propor- 
tional to the degeneracy of the medium. From high tem- 
perature phase to low temperature phase, over a nar- 
row temperature range, degeneracy of the medium drops 
rapidly, by factor of 4-5. 

Can we construct an experimental observable equiv- 
alent to s/T 3 , i.e. degeneracy of the medium? Vari- 
ation of the observable with collision energy then can 
answer, whether or not threshold energy exists for the 
confinement-deconfinement transition. We argue that 



the experimental observable, 



ratio of the d> me- 



son multiplicity ({N'*')) and cube of the <f> mesons mean 
transverse momentum ((p T )) is effectively proportional 
to s/T 3 , i.e. the de generacy of the medium. The ar- 
gument is based on three assumptions: (i) viscous ef- 
fects are small in relativistic energy heavy ion collisions, 
(ii) <j) mesons multiplicity is proportional to initial en- 
tropy density and (iii) <fi mesons mean px is proportional 
to initial temperature. Let us examine the assumptions 
in detail. At RHIC energy collisions, experimental data 
at low pt < l.bGeV are consistent with ideal hydrody- 
namics [2l| . Assumption (i) is approximately valid in 
RHIC energy. However, viscous effects can be substan- 
tial at lower SPS energy [22|, [HI • Recently Petersen and 
Bleicher [24j studied elliptic flow at SPS energy. It was 
shown that initial conditions can have substantial effect 
on development of the elliptic flow. With proper ini- 
tial conditions and with gradual freeze-out, elliptic flow 
data at SPS energy can be explained in an ideal hydro- 
dynamic model. In [25], <fi mesons transverse momentum 



spectra, in the energy range -^5=6-200 GeV, were anal- 
ysed. Ideal hydrodynamics reasonably well explained the 
spectra. Assumption of small viscous effect in relativis- 
tic heavy ion collisions seems to be reasonable at SPS 
energy also. If viscous effects are small, then assumption 
(ii) i.e. 4> meson multiplicity is proportional to initial en- 
tropy density, is also reasonable [2q . It is generally be- 
lieved that total multiplicity is proportional to final state 
entropy density. If viscous effects are small, entropy is 
not generated and initial and final state entropy is same. 
Experimental data on production cross-section of various 
particles and anti-particles show surprisingly good agree- 
ment with thermal abundances of a hadronic resonance 
gas [27|-|29(. The assumption is valid in a thermal model. 
However, some of the thermal models use 'strangeness 
under saturation (75 ) ' factor to explain the data [Til . [l5j . 
7s « 0.6 in low AGS energy. It increases with energy and 
at RHIC energy 75 w 1. The assumption may not hold 
in thermal models with strangeness under saturation fac- 
tor. The assumption (iii), <fr mesons mean px is propor- 
tional to the initial temperature is approximately valid 
in an ideal hydrodynamic model. In Fig|5J (ideal) hydro- 
dynamic model predictions for <f> mesons mean px in a 
central (0-10%) Au+Au collisions, are shown as a func- 
tion of 'spatially' averaged initial temperature. At the 
initial time Ti=0.6 fm, initial energy density is assumed 
to be distributed as [21( 



e(h,x,y) = Ei[0.75N part (b,x,y) + 0.25N coU (b, x, y)}. 

(1) 

In EqflJ b is the impact parameter of the collision and 
N part and N co u are the transverse profile of the average 
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FIG. 4: The filled circles are spatially averaged initial temper- 
ature of the fluid, obtained from an ideal hydrodynamic model 
analysis of (f> meson pr spectra in the energy range -^5=6-200 
GeV. The solid line is a fit to the average temperature. 
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FIG. 5: The experimental ratio R is plotted as a function of 
temperature. The solid line is a step function fit to the ratio. 
For comparison, we have also shown (rescaled) s/T 3 in lattice 
simulation by cheng et al. @], the dotted line. 



number of participants and average number of binary col- 
lisions respectively, calculated in a Glauber model. Initial 
fluid velocity is zero, v x (x,y) = v y (x,y) = 0. Hydrody- 
namic equations are solved with the code AZHYDRO- 
KOLKATA, detail of which can be found in We 
have used an equation of state, with a cross-over tran- 
sition from QGP phase to hadronic phase at tempera- 
ture T co =196 MeV [32]. For a set of central energy den- 
sity, £j, <fi mesons mean pr's are calculated at the freeze- 
out temperature Tp=150 MeV. In ideal hydrodynamics, 
a linear relation between (pj) and average initial tem- 
perature (Ti) is accurately observed. It is interesting to 
note that the relation (ppj oc (T/) will not be valid, say 
for the pions or for the strange mesons K + /K~. Pions 
and kaons are largely affected by resonance decay, which 
spoils the relation. Resonances do not contribute to cj) 
meson production making the relation works. 

NA49 and STAR collaboration have tabulated 4> mul- 
tiplicity and mean pr in Pb+Pb and Au+Au/Cu+Cu 
collisions in the energy range ^fs— 6.3-200 GeV [17T - I201 ] . 
In Fig(3l the experimental ratio (N^) / (p^) 3 is shown as 
a function of the collision energy. In Au+Au/Cu+Cu 
collisions, </> mesons are measured in a number of col- 
lisions centrality. Presently, we chose the most central 
ones. To account for the differences in system size, col- 
lision centrality and rapidity gap (AY) in NA49 and 
STAR data sets, 4> multiplicity is normalised by the fac- 
tor 0.5N part AY . Uncertainty in the experimental ratio 

(N^)/(pp) 3 is large due to the large experimental er- 
ror in determination of c/> meson multiplicity and mean 
Pt- For example, in low SPS energy, uncertainty in (f> 
multiplicity is as large as ~20-40%. At RHIC energy 
also, (N't') is determined only within ~ 10-20% accuracy. 
Mean px is determined more accurately, in the energy 



range -^=6.3-200 GeV, uncertainty in < pj, > is less 
than 10%. Though error bars are large, the trend of the 
ratio, R = (N*) /[.bN part AY(p^) 3 } as a function of the 
collision energy is evident. Mimicking the temperature 
dependence of s/T 3 , the ratio sharply rises from low SPS 
energy to RHIC energy. In Fig|3l the solid line is a fit to 
the ratio with an analytical form for the step function, 

R = a[l + tanh ^~/* th }. (2) 

Analytical form Eq|3J well explain the data, \ 2 /N = 
0.3. Fitted values are a = 0.017 ± 0.009, ^fs th = 
15.74±8.10GeVand A^~s th = 14.52±14.93 GeV. Thresh- 
old energy can be determined only within ~50% accu- 
racy, the width of the transition is uncertain by ~100%. 
One note that presently, no data exist in the energy range 
17.3-62 GeV, between the top of SPS energy and bottom 
of the RHIC energy. STAR energy scan programme will 
fill up the gap. Threshold energy and width of the tran- 
sition can be determined more accurately. 

One may argue that fitting the ratio (N^) /(p^) 3 by an 
analytical form for the step function is rather arbitrary, 
the ratio could as well be fitted by another form, with- 
out any threshold energy. As an example, in FigOJ we 
have shown a fit (the dashed line) to the ratio R, by a 
power law, R — A*Js B . Power law also explains the data, 
but with increased "X 2 /N = 0.92. Since the step function 
fit as well as the power law fit give x 2 /iV < 1, from 
the \ 2 analysis point of view, both the fits are equiv- 
alent and it can not be claimed that the experimental 
ratio (N^) / (plpj 3 exhibit step function like behavior. In- 
deed, one does wonder, whether it can be verified that 
the ratio (N^) / (p^) 3 do indeed corresponds to s/T 3 , i.e. 
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FIG. 6: The solid line is the square speed of sound (cj) in lat- 
tice simulation Q| . The dashed line is c? s in the model assum- 
ing the step function fit to the experimental ratio (iV^}/(py} 3 
is proportional to s/T 3 , the entropy density over cube of the 
temperature. In the inset, trace anomaly in lattice sim- 
ulation and in the model are plotted. Model trace anomaly is 
normalised by the factor K — 460. 

the degrees of freedom of the medium, as argued here. 
Note that s/T 3 is a function of temperature, not of colli- 
sion energy. Can we convert the collision energy depen- 
dence of (N^) I ' (y%) 3 to temperature dependence? As 
noted earlier, in [25| , (f> meson pt spectra in the energy 
range 6-200 GeV were analysed in an ideal hydrodynamic 
model. Ideal hydrodynamics reasonably well explain the 
Pt spectra of <j> mesons. In FigJU spatially averaged ini- 
tial temperature of the fluid as obtained in the hydro- 
dynamic analysis (25|, are shown as a function of the 
collision energy. Spatially averaged initial temperature 
depend logarithmically on the collision energy, 

<Ti>= A + Blogy/7s (3) 

with A = 0.126, B = 0.063. We use Eq|3]to convert the 
collisional energy dependence of the ratio, (N^)/(p T ) 3 
in to temperature dependence. FigJS] shows the ratio as 
a function of temperature. As a function of tempera- 
ture also, the ratio show a jump from low temperature 
to the high temperature. The ratio is still fitted with a 
step function, with T th = 202 ± 23MeV, AT th = 33 ± 37 
MeV. The threshold temperature, 1^=202 MeV, is very 
close to the cross-over temperature T co =197± 3 MeV, 
obtained in the recent lattice simulation The dotted 
line in FigJSJ is a parametric representation for the en- 
tropy density in lattice simulation [H . It has been scaled 
by a factor K = 460. Experimental ratio (AT*) / '(p T ) 3 
closely agree with lattice simulated s/T 3 and also with 
the step function fit to the ratio. Only in the temperature 
range 210-240 MeV, where data are sparse, rescaled s/T 3 



overestimate the step function fit by 10-20%. Close sim- 
ilarity between lattice simulated s/T 3 and experimental 
ratio (N^) / (p T ) 3 do support our argument that the ratio 
(N^) / (p T ) 3 is proportional to s/T 3 or the degeneracy of 
the medium. To confirm that the ratio (N^) / (p T ) 3 is 
indeed proportional to s/T 3 , using the usual thermody- 
namic relations, 

p(T) = f s(T')dT' (4) 
Jo 

e(T) = Ts-p, (5) 

we obtained pressure and energy density for the step 
function fit to the ratio (N^) / (p T ) 3 and computed square 
speed of sound ((% = ■£)• In Fig® the dashed line is c 2 , 
obtained from the step function fit to the experimental 
ratio (N^) / (p T ) 3 . The solid line is c 2 from a parametric 
representation of the lattice simulations for s/T 3 . For 
T >200 MeV, c 2 from the step function fit to the experi- 
mental ratio closely agree with c 2 in lattice simulation. c 2 s 
in lattice simulation show a dip at T ss200 MeV. c 2 from 
the step function fit to the ratio show a dip at T «180 
MeV. However, at low temperature, c 2 is comparatively 
larger in lattice simulation than in the step function fit to 
the ratio. We note that experimental data for the ratio 
{N' t ')/{p T ) 3 do not exit below temperature T «170 MeV 
and the step function fit may not be accurate at low tem- 
perature. In the inset of FigJBl the trace anomaly £ T 3p , 
is shown as a function of the temperature. The solid line 
is the trace anomaly in lattice simulation, the dashed line 
is that from the step function fit to (N^) / (p T ) 3 . Trace 
anomaly from the step function is normalised by a factor 
K = 460. Normalised trace anomaly closely agree with 
the lattice simulated value. We may note here that if the 
ratio (N^)/ (p T ) 3 is fitted with a power law, c 2 and trace 
anomaly from the fit do not resemble lattice simulation 
results. For example, square speed of sound does not 
show a dip or the trace anomaly does not show a peak. 

To conclude, we have argued that the ratio of <j> mesons 
multiplicity over cube of the mean px corresponds to 
the degeneracy of the medium produced in relativis- 
tic heavy ion collisions and can signal the confinement- 
deconfinement phase transition. The ratio constructed 
from the existing data p7l - l20j in the energy range 
\/i=6.3-200 GeV, even though have large error bars, do 
show rapid rise from low SPS energy to high RHIC en- 
ergy. Using a result of hydrodynamical analysis of <p 
mesons px spectra [25|, that the average initial temper- 
ature depend logarithmically on the collision energy, the 
collision energy dependence of the ratio is converted in to 
temperature dependence. The ratio as a function of the 
temperature closely corresponds to lattice simulations for 
s/T 3 , the entropy density over cube of the temperature. 
From a step function fit to the ratio we also have ex- 
tracted the threshold energy, \/s th — 15.74 ± 8.10 GeV, 
for the confinement-deconfinement phase transition. 
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